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Abstract 

We discuss the statistical properties of the volume of the nodal 
set of wave function for two paradigmatic model systems which we 
consider in arbitrary dimension s > 2: the cuboid as a paradigm for a 
regular shape with separable wave functions, planar random waves as 
an established model for chaotic wave functions in irregular shapes. 
We give explicit results for the mean and variance of the nodal vol- 
ume in arbitrary dimension, and for their limiting distribution. For 
the mean nodal volume we calculate the effect of the boundary of the 
cuboid where Dirichlet boundary conditions reduce the nodal volume 
compared to the bulk. Boundary effects for chaotic wave functions 
are calculated using random waves which satisfy a Dirichlet boundary 
condition on a hyperplane. 

We put forward several conjectures what properties of cuboids gen- 
eralise to general regular shapes with separable wave functions and 
what properties of random waves can be expected for general irregu- 
lar shapes. These universal features clearly distinct between the two 
cases. 



* sven. gnutzmann@nottinghani. ac . uk 



1 



1 Introduction 



We consider real square-integrable eigenfunctions $(q) of the free stationary 
Schrodinger (or Helmholtz) equation 



on an s- dimensional smooth connected compact Riemannian manifold M. 
with local coordinates q = (g^, . . . , g*) . Here, A^v; is the Laplace-Beltrami 
operator on Ai and E is an energy eigenvalue. We have set the value of the 
physical constant ^ of Planck's constant squared over twice the mass of the 
particle equal to one by appropriate choice of units. 
If Ai has a boundary we will impose Dirichlet boundary conditions. 
Compactness ensures a discrete and non-negative energy spectrum that we 
arrange in ascending order as < -Ei < £'2 < ■ ■ ■ < Ej^ < -Etv+i < ■ ■ ■ • The 
eigenfunction with eigenvalue E^ will be denoted by $7v(q) (for degenerate 
eigenvalues this requires a choice of eigenbasis and order). 
For a given eigenfunction $Ar(q) the nodal set 



consists of all interior points on the manifold where the eigenfunction van- 
ishes. For the real wave functions that are under consideration here the nodal 
set is a collection of hyper- surfaces. 

More than 200 years ago Chladni |1] visualised the vibration modes of 
plates with sand that accumulates along nodal lines. He analysed in detail the 
geometric patterns formed by the nodal lines. About 30 years later Sturm's 
oscillation theorem [2] that the n-th eigenfunction of a Sturm-Liouville dif- 
ferential operator has n — 1 nodal points may have been the first rigorous 
mathematical result concerning the nodal set of wave functions. Since this 
time the nodal set for wave functions of various types has attracted the at- 
tention of many mathematicians and physicists - and many seminal results 
followed. 

In the present work we will focus on the size of the nodal set which we 
measure by its hyper-surface volume. We denote the hyper-surface volume 
of the nodal set of the A^-th eigenfunction by at - we will refer to it just as 
the nodal volume. 

With increasing energy the wavelength becomes smaller and so does the 
typical distance between two nearby nodal surfaces. One thus expects that 




(1) 
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typical "Htv increase with the energy E]^. Comparison of nodal volumes 
between eigenfunctions at very different energies or of eigenfunctions on dif- 
ferent manifolds is possible through the dimensionless rescaled nodal volume 



(3) 



VVE^ 



where V is the volume of the manifold Ai. Note that the rescaled nodal 
volume is only defined for positive energies E^- > 0. In manifolds without 
boundary this excludes the ground state with energy Ei = which has a 
vanishing nodal set and thus "Hi = 0. 

In the mathematical literature the nodal volume has been a central object 
and the main results are best summarised by Yau's conjecture [3] that states 
that the rescaled nodal volumes are bounded from below and above for all 
smooth compact manifolds. That is 



for all > 2 where the constants ci and C2 only depend on the manifold and 
the metric. For real analytic manifolds this is a classic theorem by Donnelly 
and Fefferman [1]. For the smooth case lower bounds have been established 
[5|, [6l [7]. The general proof of Yau's conjecture remains a central problem 
in spectral theory - we refer to the recent survey by Zelditch for a more 
complete overview and additional results [S]. 

In this work we will consider how the rescaled nodal volume for a given 
manifold is distributed statistically. Such a statistical approach is well es- 
tablished for the number of nodal domains [9], [10] which revealed that shapes 
with a chaotic ray dynamics have a universal distribution which can clearly 
be distinguished from distributions for shapes with a separable Laplacian 
and thus integrable ray dynamics. For separable shapes these distributions 
can often be calculated exphcitly [9l |Tl] and share some universal features. 
Here we will consider nodal volumes for two models: 

i. The cuboid. This is a paradigm of a regular shape with separable 
Laplacian and integrable ray dynamics; 

ii. The boundary adapted planar random wave model [12]. In two dimen- 
sions statistical properties of nodal volumes and nodal densities have 
been discussed in detail for various types of random superposition of 



Cl < CTat < C2 



(4) 



3 



eigenstates pSl H [ISl [151 CZl [ISl US] • The the boundary-adapted ran- 
dom wave model is an extension of the standard Gaussian random wave 
model introduced by Berry [20] who conjectured that eigenfunction 
statistics for chaotic systems follows the predictions of the Gaussian 
random wave model in a semiclassical limit. The boundary adapted 
random wave model is able to predict systematic corrections near the 
boundary - for the two-dimensional case the effect of a boundary on 
nodal densities has been discussed in detail [121 1131 El] • We will add 
some results for higher dimensions and discuss implications for mani- 
folds with chaotic ray dynamics. 

In section[2]we will give a full derivation of limiting distributions for the nodal 
volume based on Poisson summation. Moreover we will consider boundary 
corrections to the mean of the nodal volume. In section [3] we will summarise 
Berry's results on the nodal volumes in boundary adapted random waves and 
derive some extensions. The main result of this chapter is formulated as a 
conjecture on the finite energy correction to the mean nodal volume in an 
irregular shape for arbitrary dimension. Implications of our findings to more 
general regular and irregular shapes will be discussed further in section |H 

2 The regular case: nodal volume statistics 
for an s-dimensional cuboid 

Let ai {£ = 1, . . . , = s) he the side lengths of an s-dimensional cuboid with 
volume V = Y[e=i ^i- Separation of variables leads to a unique basis of 
(normalised) eigenf unctions 



which are labelled by s positive integers ni {£ = 1, . . . , s). The corresponding 
energies are 






(6) 



and the rescaled nodal volumes are given by 
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Our first goal will be to obtain the asymptotic mean value of in a spectral 
interval € [E,E + AE] of width AE near the energy E 

{(^n)[E,E+AE] = 77 '^n X[E,E+AE]{En) ■ (8) 

NlE,E+AE] 

Here N[e,e+ae] is the number of eigenfunctions with energies in the interval 
En G [E,E + AE] and X[e,e+ae]{x) is the characteristic function of this 
interval. 

We will be interested in the asymptotic behaviour as — > oo such that 
N[e,e+ae] — )■ cxo at the same time. The choice AE = gE^^^ for some constant 
g > satisfies the above requirement and is sufficiently small such that 
any systematic change over the interval only leads to corrections which are 
much smaller than the ones we will explicitly calculate. The standard tool 
for extracting asymptotic behaviour in this setting is Poisson's summation 
formula. We will apply it in the form 

Yl /(^) = / f(^)dn - -/(O) + 2 ^ / f{n) cos{27rMn)dn (9) 

n=l ^ M=l-^0 

which is valid for sufficiently well-behaved functions f{n) (such that all sums 
and integrals converge). Poisson summation allows us to find the asymptotic 
behaviour of N[e,e+ae] = N{E + AE) - N{E) where 

NiE) = J2^(^- ^n) (10) 

=iVweyl(^)+iVosc(^) (11) 

is the spectral counting function. In the second line we have written N(E) 
as a sum of a smooth monotonically increasing function A\veyi(-E) and an 
oscillating function Nosc{E). Application of the Poisson summation formula 
reveals that N{E) is asymptotically dominated by the smooth part 

^Wo„(£) = gs"^ - ^It^B'-'"^ H- 0(£<-"^) (12) 



which is known as Weyl's law. Here = Y'(^s_^_-^'j is the volume of the s- 

dimensional unit ball and S = 2V T^f , — is the s — 1 dimensional volume 
of the surface of the cuboid. 
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Figure 1: Mean (main panels) and variance (insets) of the rescaled nodal vol- 
ume as a function of energy. Data points are averages over ^ 10^ eigenstates 
in an interval around the given energy. 



The first term in (fT2l) gives the leading growth of the number of states with 
increasing energy. The second term is the leading correction - a boundary 
effect as can be seen from the appearance of the surface volume. Each oscil- 
lating contribution is of order 0{E^^~^^^^) and thus asymptotically smaller 
than the boundary correction. 
Weyls law implies the estimate 

which includes the effect of the boundary in the leading correction. Here 
Ai? = gE^/'^ has been used to give the order of further corrections which will 
be neglected in the sequel. 

We may now apply Poisson summation to find the mean value of the rescaled 
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nodal volume 



(14) 



where 



^ _ (s - l)7r^,„2 , TT^s js - l)7r^,„i 

The three terms in the above expression for /3s have different origins in the 
asymptotic expansion of the Poisson sum. The nodal volume contains 
terms proportional to rii — 1. The terms proportional to rii dominate give 
the leading term in f lT^ and the first term of Pg. The unit shift leads to the 
second term in and the third term comes from the boundary correction to 
the spectral counting function. Note that /J^ is a positive constant for any s 
and that it decays as the dimension grows. 

Analogously one may calculate the variance which we only give to leading 
order 

Var(a„)[s,£;+AE] = ^ + ^ + ^(^ 0- (16) 

Higher moments can be obtained in a similar way. Alternatively one may 
just consider the limiting distribution 

Ps{a) = lim {6{a - crn))[E,E+AE] (17) 
for which one may derive the formal expression 

P,W^-L/d.-.n*(.-S=^) (18) 

where the integral is over the a unit s — 1-dimensional sphere and e = 
(ei, . . . ,es) is a point on the sphere (i.e. = 1). It is straight forward to 



see that these limiting distributions vanish outside the interval a G 
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P{a) 
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s = 2 
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Figure 2: Nodal volume distributions for an s-dimensional cuboid. The red 
curves show histograms obtained from 10^ highly excited eigenstates. For 
s = 2 and s = 3 the upper panels also show the limiting distributions f lT9|) 
and fl2U]) as a dashed line. 



For low dimensions (fT8|) can be given by direct integration 

4 for a G [I/tt, V2/7r], 




else; 
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V3 



^ + arctan ^ 



V6-37r^cr^ 



+ 



arctan 
arctan 



7ro-+3v'2-7r^o-^ 



47r 
V3 



else. 



(19) 



for (7 G [I/tt, V2/7r], (20) 
for (T e [v^/tt, VS/tt], 



For larger dimensions direct integration of the expression 0181) may be per- 
formed with increasing effort. In Figure [2] histograms of these distributions 
for a finite energy interval are shown that illustrate how the distribution 
changes with increasing dimension. 

Let us conclude this section with a few observations about this example which 
may be generalised to more general separable cases. Apart from points of 
high symmetry separable wave functions all have a similar local checker board 
structure. One may then expect that the qualitative behaviour of the mean 
nodal volume (including the boundary effect), the variance and the limiting 
distribution of nodal volumes will be very similar. We conjecture that the 
following features of the limiting distribution are universal 

i. The limiting distribution has compact support cr e [(Tmin? c"max] with 
Cmin > (this is consistent with Yau's conjecture). 

ii. Near a^in the behaviour is P{cr) = Ci(o" — cTramY"^ ^(o" — c"mm)- 

iii. Near (Tmax the behaviour is P{(t) = C2(o"max — cr)('^~-^)/2 ^(cTmax — cr) 
(see [TT] for an analogous singularity in limiting distributions of nodal 
counts). 

3 The chaotic case: nodal volume statistics 
of boundary-adapted random waves 

Let us now consider nodal volume statistics for s-dimensional random waves 
as proposed by Berry [20]. These are a model for wave functions in chaotic 
billiards. In order to account for boundary corrections to the mean nodal 
volume we use boundary adapted random waves following closely the anal- 
ysis presented by Berry in [12], where the mean nodal volume for random 
waves in s = 2 dimensions which satisfy a Dirichlet boundary condition along 
an infinite line was calculated (see also [131 El ESj)- the mathematical 
literature similar approaches have been used to study statistical properties 
of nodal volumes on tori [THl [T7] and spheres [THl [T^. In these cases random 
wave models have been constructed in terms of random superpositions of de- 
generate eigenfunctions and rigorous results on expected nodal volumes and 
on the fluctuations were obtained. The latter are largely consistent with the 
results obtained by Berry for planar random waves for appropriate choices 
of eigenspaces. 
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Let us now construct a Gaussian random wave model whose realisations 
are solutions of the s-dimensional free stationary Schrodinger equation ([T]) 
on the Euchdean space at energy E = k"^ with a Dirichlet condition 
$ = on the hyperplane Xg = (where x = {xi, . . . , Xg) € are Cartesian 
coordinates). We will consider solutions in the half space Xs > and thus 
refer to the plane = as the boundary. It will be convenient to use 
'dimensionless' (rescaled) coordinates, so we define R = (i?i, . . . , Rg), where 
f^'i — k X 2 • 

Let us first define the standard random waves (without any boundary 
conditions) by 

«(R) = ReJ-5^e^^-^+^<^^ (21) 
i=i 

where iij are uniformly distributed on a unit (s — l)-sphere and the phases 
(pj are uniformly distributed on [0, 27r). The Gaussian random wave model is 



achieved in the formal limit N ^ oo. The random waves (1211) do not obey any 
boundary conditions. Dirichlet boundary conditions at the boundary Rg = 
can be implemented in a straight forward way by anti-symmetrisation with 
respect to the boundary 

$(R) = i=(^R)-«(R)) (22) 

where R = R — 2Rses where = (0, . . . , 0, 1) is the unit vector in direction 
Xs- It is expected that the effect of the boundary at = becomes weaker 
when Xs becomes larger. We have chosen the normalization constant in (1221) 
so that ($^) — )■ 1 as -Rs ~^ oo where (■) refers to the average over random 
waves. 

For a given region G in the half space Xs > {Rs > 0) we may now write 
the expected rescaled nodal volume as 



kf^d 



X 



/,^(5($(R)) iWnH md'R (23) 



J,aP{R)d^R 
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where kG = {/cx : x G G} is the rescaled region and we have introduced the 
(expected) nodal density 



p(R) = (5(<l>(R)) |Vr$(R)|) . (24) 

Note that the random wave model is translation invariant with respect to 
translations parallel to the boundary Rg = such that the nodal density 
only depend on the distance Rg from the boundary 

p(R) = p{Rs) . (25) 

As — 7- oo any probability distribution the random wave model becomes 
a Gaussian process such that any probability distribution involving $ and 
its derivatives 9j$ = ^ at a point R is a multivariate Gaussian. For 
the present purpose we need P($, . . . , (9^$), and more specifically in 
P(<l> = 0,9i<l>,...,9,$). 

The relevant variances and cross-correlations can be calculated from the 
known two-point correlator of the standard random wave model which is 
given by 

s-2 /s\ J£^(Ri~R2) 

(«(R0«(R2)) = 2^ r(-) ^ ^J' . (26) 

|Ri — R2I 2 

The lengthy but straight forward calculation of the nodal density of the 
boundary-adapted random wave model can be performed by generalising 
Berry's 2-dimensional calculation [12]. We refer to [21] for details of the 
calculation which leads to the nodal density 



p{Rs) = Pbuik^ ^ F (-1 , 1 ; ^ mJ (27) 
where 

Pbuik = {6 (n(R)) |VRn(R)|) = ^"^^ (28) 

is the constant nodal density of the standard random wave model with- 
out boundary (below we will show that p{Rs) — )■ Pbuik as Rg — )■ 00) and 
F[a,b] c; x) = 2Fi(a, 6;c;x) is a hypergeometric function [22]. The nodal 
density fl27|) depends on the distance Rg from the boundary via the covari- 
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ance functions 



B ^($(R)^) =1 - r( -) (29) 



2 



1 rfi)j£(2i?,) 

Dn, =(9i$(R)2) =- - 4r^^^ (30) 
=(9.$(R)^) =- + rm^-^^ ; ^ (31) 



2i?J 

^($(R)9.$(R)) =r(-) (32) 
where J^ix) is the n-th Bessel function and the abbreviation 



M^l-?^i^. (33) 

Let us now consider the asymptotic behaviour of the nodal density f l27p as 
Rs oo. Using the known asymptotic behaviour of Bessel functions and for 
the hypergeometric function one then obtains p{Rs) Pbuik with the leading 
order smooth and oscillatory corrections given by (see [21] for details of the 
calculation) 

^ = 1 + cri?;(-i) + c°-^^!i^^i_z5ill + o(i?;(^-^)) (34) 

Pbulk Jl 2 



where 



2^+20f(, + 2)r(^) ^ ' 

r (^) 

(josc =^lL . (36) 



The oscillatory part in ( l34l) decays much slower than the smooth corrections. 
For nodal volumes one needs to integrate and we will see that the smooth 
correction will dominate over the oscillatory part. Note that f p^ has oscil- 
latory terms of order 0{R'^^^~^^) which are formally of the same order as the 
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smooth correction. 

On the boundary Rg = one finds 



p(o) iv^(s-i)r(f)^ 



< 1 (37) 



which is consistent with the expectation that a Dirichlet boundary condition 
will lead to a suppression of the nodal density near the boundary. 

We can now come back to the estimate of the rescaled nodal volume inside 
a given bounded region. For the standard random wave model the rescaled 
nodal volume is equal to the constant nodal density (ac) = Pbuik- For the 
boundary adapted random wave model there are corrections which become 
stronger close to the boundary. To be specific let G be a cylindrical region in 

of height a such that the bottom is a connected bounded s — 1 -dimensional 
region in the hyperplane Xg = 0. The volume of G is V = aS where S is the 
s — 1-dimensional hypervolume of the bottom. The rescaled nodal volume is 

^ sC,iR.),R. ^ / ^ S . (38) 

Vfc V Jo Pbuik / 

In the high energy limit E = k"^ ^ oo the rescaled nodal density converges 
to Pbuik- Indeed, p{Rs) is a bounded function that converges to pbuik as 
Rs ^ oo such that the integral j^"' ^^^p|,~^''""' d-^s is at most of order o{k) 
and the correction term is at most o(l) for k — )■ oo. The form of the leading 
order correction can be obtained from the asymptotic expansion of the nodal 
density (134|) . In the 2-dimensional case one finds 

r P^^^l^P^dR, ~ - /2 + Cr r ^2'dR2 (39) 

Jo Pbuik Jl 

~-/2 + Criog(A;a) (40) 

where h = - - 1 - Cf^0{R2 - l)i?2 ^) dR2 is a constant term. 

For any higher dimension s > 3 we may define the constant 

'piRs) 



Pbuik 



l]dRs . (41) 



Altogether we obtain 



Pbuik (1 - + for ^ = 2, .421 

^ " ^ Pbuik (1 - f ^ + 0(1/P/^)) . ^''^ 
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By numerical integration one obtains the coefficients ~ 0.758 and I4 = 
0.645. 

Let us also mention that Berry has shown that (cr|.) — (cg)^ = (^^^) 
[T2] in s = 2 dimensions which implies that the distribution of nodal volumes 
P(a) for a finite energy interval is very narrow and converges to 

P(a) = 5(a-pbuik) (43) 

as E = k'^ ^ 00. The same scaling also applies to a random wave model on 
the 2-sphere as shown by Wigman in [19]. The latter work by Wigman also 
states that for s > 3 one may show that the rescaled nodal volume variance 
for the s-sphere is bounded by 0(k~^^~^^). This implies that fl43l) applies to 
all dimensions for random waves on spheres. As curvature effects should not 
change these scalings one expects the limiting distribution P3l) also in the 
case of the Euclidean random waves with boundary considered here. 

At the end of this section let us come back to more general manifolds with 
chaotic ray dynamics. Berry has introduced the standard Gaussian random- 
wave model without boundary as a model for wave functions in this case 
and he also proposed to include boundary effects with the boundary-adapted 
random-wave model. In [12] he conjectured that for s = 2 the corrections 
of the boundary-adapted random wave model should also apply to chaotic 
billiards. Here we extend his conjecture to arbitrary dimensions, i.e. the 
asymptotic behaviour of the nodal volume is described by ( H2l) including the 
leading order correction terms and no free parameters as V should be replaced 
by the volume of the manifold Ai and S is the hypervolume of its boundary. 
For higher order corrections one may need to include curvature effects which 
are not represented in the random wave model. In the case s = 2 terms of 
order O (k^^) in the nodal volume of random waves in a cylindrical region 
explicitly contain the height a which is not well defined when one tries to 
translate the result to general billiards |12] . 

4 Signatures of wave chaos and integrability 

In the previous two chapters we have given a detailed account of nodal vol- 
ume statistics for two paradigmatic systems in arbitrary dimension. The 
s-dimensional cuboid is a paradigm for a regular shape with separable wave 
functions while planar random waves are a model for wave functions on irreg- 
ular shapes with chaotic ray dynamics. In both cases we have found limiting 
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distributions for the rescaled nodal volume. The limiting distributions in the 
two paradigms have a very different character: for random waves we find a 
delta-function while the cuboid's limiting distribution has a finite support. 
It is consistent with Berry's random wave conjecture that the distribution 
for all irregular shapes will be a delta function at a = pbuik (see fl28|) ) which 
only depends on the dimension and nothing else. For regular shapes other 
than the cuboid we may refer to the analogous calculations for nodal count 
distributions [9l [11] which showed that many features such as the types of 
singularities near the upper and lower end of the support are universal. We 
conjecture that the same is true for the nodal volume distributions studied 
here though the actual values for the upper and lower end of the support may 
be system dependent. One clear signature of irregular versus regular shapes 
is the behaviour of the fluctuations of rescaled nodal volumes as the energy 
increases. For irregular shapes one expects that the variance decreases while 
it remains finite and bounded for regular shapes. 

Another interesting difference between regular and irregular shapes in di- 
mension 2 is the different nature of the boundary correction to expected 
rescaled nodal volumes - these go like 1/k for the cuboid (and, conjecturally, 
for other regular shapes with separable wave functions) and like log k/k for 
the random wave model. It has been shown by numerical computations that 
the boundary adapted random wave model gives an accurate account of the 
nodal density near a boundary [T3] and it was conjectured by Berry that this 
may be seen in any chaotic billiard. 

Let us also note, that the average rescaled nodal volume for the cuboid is 
always larger than the one for a random wave. This is expected as the 
nodal surfaces of a separable function generally intersect in a checker-board 
structure while nodal intersections are avoided in random waves [22] ^ this 
effective repulsion of nodal surfaces leads to the decrease in the expected 
nodal volume. Again, the same decrease can be expected for eigenf unctions 
in irregular shapes when compared to separable eigenfunctions of a regular 
shape. 

References 

[1] E.F.F Chladni, Die Akustik (Breitkopf and Hartel, Leipzig, 1802). 

[2] J.C.F. Sturm, Sur les equations differentielles lineaires du second or- 



15 



dre, J. Math. Pures et Appl. de Liouville 1, 106 (1836); Sur une classe 
d' equations differentielles partielles, J. Math. Pures et Appl. de Liouville 
1, 375 (1836). 

[3] S.T. Yau, Survey on partial differential equations in differential geome- 
try. In: Seminar on Differential Geometry, Ann. of Math. Stud. 102, 3 
(Princeton University Press, 1982). 

[4] H. Donnelly, C. Fefferman, Nodal sets of eigenf unctions on Riemannian 
manifolds, Invent. Math. 93, 161 (1988). 

[5] J. Briining, D. Gromes, Uher die Ldnge der Knotenlinien schwingender 
Membranen, Math. Z. 124, 79 (1972). 

[6] J. Briining, Uber Knoten Eigenf unktionen des Laplace- Beltrami Opera- 
tors, Math. Z. 158, 15 (1978). 

[7] T.H. Colding, W.P. Minicozzi II, Commun. Math. Phys. 306, 777 
(2011). 

[8] S. Zelditch, Eigenf unctions and nodal sets. larXiv: 1205.28121 [math.SP]. 

[9] G. Blum, S. Gnutzmann, U. Smilansky, Nodal Domain Statistics: A Cri- 
terion for Quantum Chaos, Physical Review Letters, 88, 114101 (2002). 

[10] E. Bogomolny, C. Schmit, Percolation Model for Nodal Domains of 
Chaotic Wave Functions Phys. Rev. Lett. 88, 114102 (2002). 

[11] S. Gnutzmann, S. Lois, On the nodal count statistics for seperable 
systems in any dimension, accepted for publication in J. Phys. A; 
larXiv:1208.2T20l [math-ph]. 

[12] M.V. Berry, Statistics of nodal lines and points in chaotic quantum bil- 
liards: perimeter corrections, fluctuations, curvature, J. Phys. A. 35, 
3025 (2002). 

[13] M.V. Berry, H. Ishio Nodal densities of Caussian random waves satisfy- 
ing mixed boundary conditions, J. Phys. A. 35, 5961 (2002). 

[14] M.V. Berry, H. Ishio Nodal-line densities of chaotic quantum billiard 
modes satisfying mixed boundary conditions, J. Phys. A. 38, L513 (2005). 



16 



M. Dennis, Nodal densities of planar gaussian random waves, Eur. Phys. 
J. Special Topics 145, 191 (2007). 

F. Oravecz, Z. Rudnick, I. Wigman, The Leray measure of nodal sets for 
random eigenf unctions on the torus, Ann. Inst. Fourier 58, 299 (2008). 

Z. Rudnick, I. Wigman, On the volume of nodal sets for eigenf unctions 
of the Laplacian on the torus, Ann. Henri Poincare 9, 109 (2008). 

I. Wigman, On the distribution of the nodal sets of random spherical 
harmonics, J. Math. Phys. 50, 013521 (2009). 

I. Wigman, Fluctuations of the nodal length of random spherical har- 
monics Comm. Math. Phys. 298, 787 (2010); erratum Comm. Math. 
Phys. 309, 293 (2012). 

M.V. Berry, Regular and irregular semiclassical wave functions, J. Phys. 
A 10 2083 (1977). 

S. Lois, thesis to be submitted to the University of Nottingham for the 
degree of Doctor of Philosophy. 

M. Abramowitz and I. Stegun, Handbook of Mathematical Functions, 
(Dover Publications, New York, 1965) 

A.G. Monastra, U. Smilansky, S. Gnutzmann, Avoided intersections of 
nodal Unes J. Phys. A 36, 1845 (2003). 



17 



